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ABSTRACT 



We first briefiy review how we investigate the modes of oscillation trapped 
c§ ! within the inner region of accretion disks by the strong-field gravitational prop- 

erties of a black hole (or a compact, weakly-magnetized neutron star). Then we 
^ ! focus on the ' corrugation' (c)-modes, nearly incompressible perturbations of the 

inner disk. The fundamental c-modes have eigenfrequencies (ordered by radial 
mode number) which correspond to the Lense-Thirring frequency, evaluated at 
the outer trapping radius of the mode, in the slow rotation limit. This trapping 
radius is a decreasing function of the black hole angular momentum, so a signif- 
icant portion of the disk is modulated only for slowly rotating black holes. The 
eigenfrequencies are thus strongly increasing functions of black hole angular mo- 
mentum. The dependence of the eigenfrequencies on the speed of sound within 
(or the luminosity) within the disk is very weak, except for slowly rotating black 
holes. 
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Subject headings: accretion, accretion disks — black hole physics — gravitation 
— hydrodynamics — relativity 

1. Introduction 

For twenty years [beginning with Kato & Fukue (1980)], it has been known that general 
relativity can trap normal modes of oscillation near the inner edge of accretion disks around 
black holes. The strong gravitational fields that are required can also be produced by neutron 
stars that are sufficiently compact (with a soft equation of state) and weakly magnetized to 
produce a gap between the surface of the star and the innermost stable orbit of the accretion 
disk. Although we shall not explicitly consider such neutron stars here, the results obtained 
will also apply to them to first order in the dimcnsionlcss angular momentum parameter 
a = cJ/GM^, since their exterior metric is identical to that of a black hole to that order. 

These modes of oscillation provide a potentially powerful probe of both strong gravita- 
tional fields and the physics of accretion disks, since they do not exist in Newtonian gravity. 
In addition, their frequencies depend upon the angular momentum as well as the mass of the 
black hole. The fractional frequency spread of each mode depends upon the elusive viscosity 
parameter of the accretion disk. 

The subject of 'relativistic diskoseismology' has recently been reviewed by Wagoner 
(1999). In this paper we shall focus on the ' corrugation' (c) modes, previously studied (less 
generally) by Kato (1990, 1993), Ipser (1994, 1996), and Perez (1993). We shall briefly 
compare these modes with the 'gravity' (g) modes and the pressure (p) modes. 

2. Basic Assumptions and Equations 

2.1. Structure of the unperturbed accretion disk 

We take c — 1, and express all distances in units of GM/c^ and all frequencies in units of 
c^/ GM (where M is the mass of the central body) unless otherwise indicated. We employ the 
Kerr metric to study a thin accretion disk, neglecting its self-gravity. The stationary {d/dt = 
0), symmetric about the midplane z = 0, and axially symmetric {d/dip=0) equilibrium disk 
is taken to be described by the standard relativistic thin accretion disk model (Novikov & 
Thorne 1973; Page & Thorne 1974). The velocity components v'' = = 0, and the disk 
semi-thickness h{r) ~ Cg/fl -C r, where Cs{r,z) is the speed of sound. The key frequencies. 
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associated with free-particle orbits in the disk, are 

n{r) = (r3/2 + a)-i , 
Ox(r) = 0(r) (1 - 4a/r3/2 + Sa^r^Y^^ , 
K{r) = n{r){l-6/r + 8a/r^/^ -Sa^r^Y^^ ; (2-1) 

the rotational, vertical epicyclic, and radial epicyclic angular velocities, respectively. The 
angular momentum parameter a is less than unity in absolute value. 

The inner edge of the disk is at approximately the radius of the last stable free-particle 
circular orbit r = rj(a), where the epicyclic frequency /t(rj) = 0. This radius rj(a) is a 
decreasing function of a, from rj(— 1) = 9 through ri(0) = 6 to ri(l) = 1. So all the relations 
we use are for r > r^, where K{r) > 0. Note, in particular, that 

Q(r) > n±{r) > K{r) , a > ; n±{r) > Q(r) > K{r) , a < . (2-2) 

For a = 0, i. e., a non-rotating disk, fl{r) = fl±{r) > K{r). 

To simplify the analysis, we here consider mostly barotropic disks [p = p{p), vanishing 
buoyancy frequency; a generalization to a small non-zero buoyancy is described briefly in 
section 4.3]. In this case hydrostatic equilibrium provides the vertical density and pressure 
profiles 

p = po(r)(l - y'y/^^-'^ , p = Po(r)(l y2)^V(r-i) (r > 1) . (2-3) 
The disk surfaces are at y = ±1, with y related to the vertical coordinate z by 

z lr-1 

The adiabatic index F = 4/3 within any radiation pressure dominated region of the disk, and 
F = 5/3 within any gas pressure dominated region. More information on the unperturbed 
disk can be found in sections 4.3 and 5.1. 



2.2. Equations for the disk perturbations 

To investigate the eigenmodes of the disk oscillations, we apply the general relativistic 
formalism that Ipser & Lindblom (1992) developed for perturbations of purely rotating 
perfect fluids. Neglecting self-gravity, it allows one to express Eulerian perturbations of all 
physical quantities through a single function 5V oc 5p/p which satisfies a second-order partial 
differential equation. Due to the stationary and axisymmetric background, the angular 
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and time dependences are factorized out as 5V — y(r, 2;) exp[i(m0 + at)], where a is the 
eigenfrequency. Then the assumption of strong variation of modes in the radial direction 
(characteristic radial wavelength <C r) ensures the WKB separability of variables in 
the partial differential equation for the functional amplitude V{r,z) = Vr{r)Vy{r,y). The 
function Vy varies slowly with r. The resulting ordinary differential equations for the vertical 
{Vy) and radial (VJ.) eigenfunctions are [see Nowak & Wagoner (1992) and Perez et al. (1997) 
for details] : 



2y dVy , 2u;l 



(fVr 



r 

d 

dr 



1 dy 



r - 1 

dVr 



1 + 



,2/, .2 2\ 



dr 



+ a'{u;'-K:') 1- 



(2-4) 
(2-5) 



The coefficient a{r) coincides with l/cs(r, 0) up to a relativistic factor of order unity which 
varies slowly with radius. 

Together with the appropriate homogeneous boundary conditions [discussed by Perez et 
al. (1997)], these equations generate the vertical and radial eigenvalue problems, respectively. 
The radial conditions depend on the type of mode and its capture zone, as discussed in the 
following sections. The coefficient a, the vertical eigenfunction Vy and eigenvalue (separa- 
tion function) \E' vary slowly with radius, as does the dimensionless ratio of the corotation 
frequency a;(r) and fl±{r): 



a;*(r) = a;(r)/Qj_(r) , a;(r) = cr -|- mfi(r) . 
(Note that ^' and a;* depend also on the eigenfrequency a.) 



(2-6) 



Our goal is to determine the vertical eigenvalues \l/(r) and the corresponding spectrum 
of eigenfrequencies a from the two eigenvalue problems for equations (2-4) and (2-5). 



3. Classification of Modes 

One can see from the radial equation (2-5) that a mode oscillates in a domain of r > 
where the quantity (o;^ — «;^)(1 — ^/loI) is positive. The first factor here (for a in the allowed 

range) changes sign twice at points r = r±(m, a, a), Vi < r- < r^, determined by equations 
(1) and (2-6). This factor is negative between r_ and r+ and positive otherwise. In Figure 
1 radii r± are plotted as functions of a for a few values of m and a. Thus the classification 
of modes is based on the behavior of the second factor, containing the eigenvalue. 

There are essentially three possibilities, corresponding to the following types of modes: 
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• g-modes: "if > (typically capture zone r_ < r < r+. 

• p-modes: ^ < (typically <C) cu^; capture zone < r < r_ or r > r+. 

• c-modes: ^ = a;^, at least at one value of the radius. 

The capture, or trapping, zone is where the mode oscillations occur. 

The g (inertial-gravity)-modes are centered on the radius where r_ — r+, correspond- 
ing to the maximum eigenfrequency \a\. The corotating eigenfrequencies corresponding 
to the lowest mode numbers are close to K{rm)- The p (inertial-pressure)-modes involve the 
uncertain physics at either the inner or outer radius of the disk. 

The c(corrugation)-modes which have so far been investigated are of a special sort; 
namely, "if — cul ~ throughout their whole capture domain. Only such c-modes are studied 
in this paper. They are typically non-radial (m = ±1) vertically incompressible waves near 
the inner edge of the corotating disk that precess around the angular momentum of the 
black hole. Their fundamental frequency is shown below to coincide with the Lense-Thirring 
(Lcnse & Thirring 1918) frequency (produced by the dragging of incrtial frames generated 
by the angular momentum of the black hole) at their capture zone boundary rc(m, a) in 
the appropriate slow-rotation limit. This interesting fact is one reason we are concentrating 
on the c-modes in this paper. The other reason is that g- and p-modes have been rather 
extensively studied by Perez et al. (1997) and Nowak & Wagoner (1991), respectively [and 
in earlier papers cited therein and in a brief summary of results by Wagoner (1999)]. 

4. Theory of C— modes and the Lense— Thirring Frequency 

4.1. Vertical eigenvalue problem: a selection rule for the axial and vertical 

mode numbers 

According to the definition of the c-modes we are interested in, we set 

*(r)M^(r) = 1 - x{r) , |x(r)| « 1 , (4-1) 
and write the vertical equation (3) as 

with the boundary condition |\^(±1)| < oo. The last term containing the slowly varying 
function x(r) is considered a small perturbation. We then immediately see that the unper- 
turbed equation coincides with that for the Gegenbauer polynomials (Bateman & Erdely 
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1953) C^(y), A = (3 - r)/2(r - 1) > -1/2, j = 0, 1,2, . . .; which are its only solutions 
regular at y = ±1. Therefore the coefficient of Vy must coincide with the corresponding 
eigenvalues: 2ul/ (F — 1) = j{j + 2A). However, the left hand side of this equality is a (slowly 
varying) function of r, so it can be at best valid only 'to the main order'. To clarify the 
meaning of that, we make use of the definition (2-6) of cj* and the fact that the difference 
between fl and fl± is typically very small, except for values of a close to 1 and simultaneously 
r close to ri ~ 1. So (anticipating that |cr| <C fl± in the same regime) we introduce a small 
parameter e by setting 

0^*^ =m(l + e), e = e(r,a,a) = '-, |e(r)| < 1, (4-3) 

and expand the eigenvalue and eigenf unction in it: 

X = exi + 0(6^) , Vy^V' + eV' + 0{e') . (4-4) 
To zero order we therefore have 

V\y) ^ V^'iy) = C}{y) (4-5) 

and 2m7(r - 1) = j{j + 2A), or 

2m2 = f (r-l)+j(3-r) , m = 0,±l,±2,...; j = 0,l,2,.... (4-6) 

Note that j is the number of vertical eigenfunction nodes through the disk. 

Equation (4-6) is the criterion for existence of the studied c-modes. It is also a selection 
rule for the angular mode numbers m and j because, depending on the value of F, the modes 
may not exist for any m and j. On the other hand, given a rational F, there might exist 
many integer values of irP' and j satisfying equation (4-6). (Using some results from number 
theory [Hua (1982)], it is possible to show that equation (4-6) has a finite set of solutions if 
and only if 2pq = k"^ for some integer k, with coprime p, q from F = 1 +p/q. For instance, we 
have an infinite sequence of c-modes for both F = 4/3 and F = 5/3, but only a finite number 
of them if, say, F 3/2.) However, the only solutions independent of F are j = m = (which 
does not provide radial mode trapping and must thus be discarded) and 

J = 1 = m% (4-7) 

which specifies the fundamental c-mode. Its vertical eigenfunction to the main order is 
linear in the vertical coordinate, oc y. Some other possible modes are = 4, j — 3; 

= 25, j — 10 for F = 4/3; and = 16, j = 6 for F = 5/3. To be specific, from now on 
we assume m > 0, a > 0, unless indicated otherwise. 
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When equation (4-6) is satisfied, the first correction xi to the vertical eigenvalue is found 
in a standard way for the discrete spectrum perturbation, from the solvability condition of 
the problem for V^: 

XI - Xi(m,j, r) = ^^^\,'^yJyo^ > 0, (4-8) 

where (■, ■) is the proper scalar product with the Gcgenbauer polynomial weight (1 — y'^)'^ 
From equations (4-1) and (4-4), the eigenvalue is then obtained as 

^ir)/u;lir)^l-exi + Oie') . (4-9) 

The integrals in equation (4-8) can be expressed in terms of the Euler gamma function in 
the general case and are elementary to calculate for j = 1, giving 

Xi = 3r-1 (m^=j = l). (4-10) 

The vertical eigenvalue and eigenfunction for all c-modes are completely determined to the 
needed order by equations (4-4), (4-5), (4-6), (4-8), and (4-9). 



4.2. Radial behavior of the vertical eigenvalue: the Lense-Thirring frequency 
at the capture zone boundary is the eigenfrequency 

We now need to study the behavior of the vertical eigenvalue as a function of the radius. 
From our definitions of the quantities involved, it follows that for a > the product 

Q^(r)e(r) = Q(r) - n^{r) + (a/m) 

is a decreasing function of the radius, because such is the difference Q{r) — 0.±{r) [see the 
basic expressions (2-1)]. Thus in this case e(r) has at most one root (zero) between and 
^maa;(the outcr radius of the disk), since fl±{r) > 0. If e(rj) < 0, then e(r) remains negative 
in the whole range ri < r < Vmax, or, equivalently, ^(r)/a;^(r) > 1 there. However, in such 
a case we are dealing, by definition, with a g-mode rather than a c-mode. 

Hence we have to assume that e(rj) > 0, which converts into a restriction on the eigen- 
frequency, a + m[fl{ri) — fi±{ri)] > 0. Then at any value of the root rg, equation (4-3) 
gives 

(T = -m [Q(rc) - fl±irc)] , n < Tc < r_ . (4-11) 

The second inequality here is proved in the following way. The definition of r_ given in 
section 3 implies (^(r_) = ±K(r_), with the plus sign for m > and the minus sign for 
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m < 0. Hence from the definition (4-3) of the parameter e(r) we derive, using equation 
(2-2): 

.(r_) = -l + ^^ = -l + ^i^<0. (4-12) 

So e(r) is positive at r = r^, negative at r = r_, and the single root Vc of this function is 
between these two points. Thus \I/(r)/ciJ^(r) < 1 for r-j < r < Tc, and ^'(r)/a;^(r) > 1 for 
Tc < r < r_. So the c-mode is trapped in < r < Tc, and Tc is the outer boundary of the 
capture zone. 

Expression (4-11) for the c-mode eigenfrequency in the case = 1 has the same 
structure as the result obtained by the approximate dispersion relation (Kato 1990, 1993), 
with a significant advantage that the value of the radius involved in it is uniquely specified as 
the right capture zone boundary of the mode in question. Of course, now it is also rigorously 
derived from the c-mode definition. 

Moreover, for r > ri{a) > 1, \a\ < 1 is smaller than any positive power of r. In 
particular, a/r^^"^ <S 1, except for a very close to 1 (where some of the above considerations 
do not apply, anyway; e. g., e(rj) — > oo as a — > 1). Using this fact, from equations (2-1) and 
(4-11) we obtain 



2a 



—m— 

^' c 



l + O 



a 



(4-13) 



In physical units, the quantity 2a/ is the Lense-Thirring frequency 2GJ / {c^r^). Thus, 
to lowest order in a/{rc)^/'^^ the c-mode eigenfrequency \a\ is an integer multiple of the 
Lense-Thirring frequency at r = Vc- In particular, the eigenfrequency of the fundamental 
c-mode (m? = j = 1) coincides with the Lense-Thirring frequency at its outer trapping zone 
boundary. For ajrl^'^ < a/r\^'^ <S 1, the fundamental c-mode is a low frequency oscillation. 

We can now replace the eigenfrequency cr by its dependence on Tc given by equation 
(4-11), 

( \ \0(\ 0( ( \} ( \ [0(r)-0x(r)]-[i7(r,)-0^(re)] 
u{r, Tc) = m [U{r) - Q{rc) + ^±{rc)\ , e(r, Vc) = ihjr) ' 

(4-14) 

and use Tc as a spectral parameter in the radial eigenvalue problem. Its eigenvalue will 
fix unambiguously the value of both the capture zone boundary Tc and the eigenfrequency 
a. Unhke other recent suggestions regarding the manifestation of the Lense-Thirring effect 
in accretion disks around black holes (Cui, Zhang & Chen 1998), the value of the radius 
corresponding to a Lense-Thirring frequency is not a (loosely constrained) parameter, but 
is predicted by the theory. 

Note that for a < the product ^}±{r)e{r) is an increasing function of the radius, so in 
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this case we have to require that e(rj) < 0. However, as seen from equation (4-12), e(r_) is 
still negative, so the single root Vc of e(r) is greater than r_. In fact, it is also greater than r_|_, 
because of the expression similar to equation (4-12), e(r+) = —1 — K(r+)/ [|m|r2_L(r+)] < — 1. 
Hence the capture zone is r_ < r < r+, and again wc arc dealing with a g-mode rather then 
a c-mode. Therefore, there are no c-modes of the studied type for counter-rotating disks. 
Accordingly, only the case of corotation is investigated in the sequel. 

Finally, both the co- and counter-rotating c-modes may have, in principle, a second 
oscillation domain, r_ < r < r_|_ and < < r, respectively. (Note that the coefficient of 
Vr in the radial equation (2-4) is also positive there.) These domains, however, arc situated 
at rather large values of the radius (up to infinity in the second case). But expression (4-3) or 
(4-14) shows that for a given value of Tc and large r, e(r) oc r^/^ ^ 1, in severe contradiction 
with the definitive property of the modes in question. Even |e(r+)| > 1, as seen from the 
above. In effect, the c-modes of the studied type can exist only in a corotating disk and are 
trapped in its innermost region. 



4.3. Vertical eigenvalue problem for slightly buoyant disks 



The above theory can be readily extended to the case of a disk with a small buoyancy. 
It can be characterized by a modified pressure dependence [compare with equation (2-3) and 
with formula (4.16) from Perez et al. (1997)] 



p = Po(r)(l-2/Y/("-^)[l-/(r,|/)], 



(4-15) 



where the even in y and slowly varying with r perturbation / is small, \f{r,y)\ <S 1. To 
lowest order this results in the following change of the vertical equation (4-2): 



dy'^ 



_2y_ 
r- 1 



^ + ^^[^-ix-S2){l-y')]Vy = 0; (4-16) 



where 



2T 



6, = 6,if)ir,y) 



r-1 



2r 



■X 



V-fA-iV/'-./l- ' 



r -1" J \ / 

S2^S,if)ir,y)^'^lf' + ^~^ 



3- 



2?/ 



(^/'-(^ + ^)/" + 0-5(1 -yV 



Here and below the prime means a derivative in y, and it is easy to check that the above 
expressions remain finite at the midplane y — because f{r,y) is even in y. According to 
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the definition (4-1), x(r) = 1 — ^{r)/u;^{r), and it is small for the studied c-mode, so there 
are three possible situations depending on how the new small corrections compare to it. If 
the buoyancy-induced terms and 62 are much smaller than i. e., than our basic small 
parameter e from the equality (4-3), then nothing changes in the main approximation. In 
the opposite case the very existence of the modes we are looking for is questionable, and the 
whole theory should be changed essentially. Thus the only interesting possibility is when 
both perturbations are of the same order. This assumption can always be made for a given 
degree of buoyancy. In such a case by the same standard technique of discrete spectrum 
perturbations we obtain, instead of the relation (4-4), the form 

X = X{r, a, m, j, T) = Xie(r) + xi + • • • , (4-17) 

where Xi{^,j, T) is still given by expression (4-8) and m, j and F are related by the selection 
rule (4-6). We obtain 

V - .~(rn , V) - ((^ " y')W)V', - - l)(2m^)-^(c5,(/)(l - y^)[V']\ V^) 
Xi - XU',", J- ; - {{1 - y'^)V^,V^) 

(4-18) 

When calculating this, one should replace cul by in the above expression for 62- The 
unperturbed eigenfunction = Vj^{y) and the scalar product (•, •) are defined in section 4.1. 

As for the radial behavior of the vertical eigenvalue studied in section 4.2 for zero 
buoyancy, it might not change qualitatively due to the new term xi, in which case only the 
value of Tc, specified now by the equation x(rc) = instead of e(rc) = 0, changes somewhat. 
Otherwise the whole mode capture zone changes. For instance, it turns into the interval 
(rj, r_) if x(r) remains positive throughout it. 

To understand the results of such possible changes, as an example let us treat a buoyancy 
perturbation which might be considered a typical one for thin disks, 

f{r:y)=^{r)y' . (4-19) 
The buoyancy corrections then become (for — 1) 

Si = ar)^-^y, S2 = ar)^^, (4-20) 

and the resulting expressions simplify greatly. In particular, for the fundamental c-mode the 
formula (4-9) generalizes to 



*(r)M'(r) = 1 - 



(3r - . <5£Z|^,„- 



+ ..., (m2=j = l). (4-21) 



It is clear that the radial behavior depends on the value of F and on how ^(r) compares to 
e(r). We consider only non- buoyant disks in the sequel. 
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5. Radial Eigenvalue Problem: Determining the Capture Zone and the 

Eigenfrequency 

5.1. Formulation of the problem 

Recall that we are now dealing with the case a > 0, in which case the c-modes are 
trapped near the inner edge of the disk, < r < Tc- We require the modes to decay within 
Tc < r < r_; this requirement provides the boundary condition at r = Tc for the radial 
equation (2-5). The boundary condition at the inner edge is not easy to specify, so we 
employ a general linear combination. 



cos — sm 6 Vr 

dr 



, < ^ < I ; (5-1) 



keeping the lack of knowledge parameterized by ^ e [0, 7r/2]. 
We introduce a new independent variable 

r(r) = r U{r') - ^{r')] dr' , = r{r,) , (5-2) 

which increases monotonically from zero to Tc when r goes from to Tc, so the inverse 
transformation is well defined. The new variable allows us to rewrite equation (2-5) in the 
WKB form: 

d'^Vr , Xiea^ 



+ 7^2Vr = . (5-3) 



The boundary condition (5-1) in terms of r becomes: 

u^n) cos ^ ^ - sin e K = . (5-4) 

a,T T=o 

The coefficient in front of Vr in equation (5-3) is obtained from its general expression (4-9) 
to lowest order in e. For convenience, we assign a special notation to it, 

Q(r) = ^ . (5-5) 

All factors here can be considered as functions of r via the dependence r = r(T) inverse to 
the function (5-2). 

We now need to describe the behavior of a(r), which is determined by the fact that this 
function is essentially the inverse speed of sound at the midplane, a{r) oc cj^(r, 0). Since the 
pressure at the inner edge vanishes in a typical disk model [p oc (r — r^)'^, with < A; < 1] if 
the torque does, we obtain a{r) oc (r — n)'^ with fj, — (T — l)k/2V. If even a small torque 
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is applied to the disk at the inner edge, which is not unhkely in practice, the pressure and 
the speed of sound become nonzero at r = rj. Hence a{r) loses its weak singularity there. 
This situation is included by assigning /i — 0. 

It is also known that the speed of sound drops as some power at large radii. Incorporating 
all this information, we are now able to specify the structure of the function a{r) as 

a(r)^7(r) "''^^^ , < < 1/2 , u > , (5-6) 

where 7(r) is some function bounded from above and away from zero, varying slowly with 
radius, and tending at infinity to a limit 'y{oo) = 7oo- A typical value of /i corresponding to 
a simple root {k = 1) of the pressure at and the value T = 5/3 is /i = 2/5. It is also true 
that u < 1 in most disk models. 

Our radial eigenvalue problem is thus completely defined. Therefore we proceed to 
determine its eigenvalues. 



5.2. WKB solution of the radial eigenvalue problem 

We apply the WKB procedure to the eigenvalue problem specified by equations (5-3), 
(5-4), and the decay condition for r > Tc. This seems to be the only effective analytical 
approach. It is validated partly by the fact that despite e(r) being small, the coefficient 
Q{t) is still rather large everywhere except close to Tc because a{r) is large. This follows 
from equation (5-6) since 7 ~ c/cg, with c the speed of light and Cs(^ c) the speed of sound 
[for estimates of the WKB approximation accuracy see Olver (1982), for example]. Thus we 
represent the solution as 

Vr OC Q-'/^(t) cos [$(t) - $e] , ^(t) = F Q^^^t') dr' (5-7) 

Jo 

in the whole capture zone < r < Tc, except small vicinities of its boundaries t — and 
T — Tc- Here $c is a constant phase, so far unknown. 

The right boundary r = Tc proves to be a turning point of equation (5-3). As shown 
in section 4.2, e(r), and hence Q(r), change sign there, being positive on the left of and 
negative on the right of it, with e(rc) = Qir,-) = 0. Near Tc wc have Q{r) = —q^ (r — Tc) + 
with = ~Q' (tc) > 0. Hence equation (5-3) reduces to the Airy equation, whose solution 
decaying for r > Tc is given by the Airy function of the first kind. 



Vr OC Ai {qy^ (r - Te)) 



(5-8) 
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Matching in the usual way the negative infinity asymptotics of the function (5-8) with the 
T — > Tc — hmit of equation (5-7), we obtain the following expression for the phase $c (where 
n is an integer related to the number of radial nodes of the mode) : 

r Qy^(r)dT + j-Tm . (5-9) 
Jo 4 



Now, close to the left boundary r = 0, according to equations (5-5), (5-6), and (5-2), 
the coefficient Q{t) has a weak singularity, 

Q{T) = q,T-'^ + ..., g,(re,a)^Xi7'(rOr?'^+'^^Kr,)|^(^-'')e(r,)>0, (5-10) 

and equation (5-3) becomes 

Its exact solution satisfying the boundary condition (5-4) is given in terms of Bessel functions 
as 



sin^l^/^J^(Ari-^)+a;2(r,)cos^^^J_C {Xr'-^) 



(5-11) 



where for brevity we have introduced the notations 



2(1-/.) 



r± = r(i ± c) 



(5-12) 



\^{z) is the Euler gamma-function.] Using the large argument expressions of the Bessel 
functions, it is straightforward to write the 'far field' (Ar^"'^ ^ 1) asymptotics of the solution 
(5-11) as 

Vr oc T^"'^ COS (XT^-^' - ^ + , (5-13) 

where 



r_ 



^i{fc, o) = arctan ( -— tan — 1 , Z±{rc, a) = —ou (rj) ( — ) cos 6 ± sin 6 . (5-14) 



2C 



The last step is to match the asymptotic expression (5-13) with the r — > -|-0 hmit of 
the WKB solution (5-7). In view of the expression (5-10), the latter is easily found to be 

Vr OC COS {Xr^-f" - $c) , 

so the match requires that 
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Employing the expression (5-9), we then arrive at the desired eigenvalue equation 

f Q^'^{T)dT ^nn . (5-15) 
Jo 

For a given set of all the parameters, including a and n, any solution of this equation 
specifies the corresponding c-mode boundary according to the relations (4-14) and (5-2). 
The eigenfrequency cr of the mode is then found from Tc by the general formula (4-11). So 
we now study the roots Tc of equation (5-15). 



5.3. The c-mode eigenfrequency spectrum is at most finite 

It is now convenient to work directly with rather than with Tc- So first of all we 
rewrite the eigenfrequency equation (5-15) in the following form: 

r{r^, a) = $(re, a) + I{r^, a) = n , n = 0, 1, 2, . . . ; (5-16) 

where ^ 

$(r„a) = -$,(re,a) , (5-17) 

TT 

7(re, a) = i ^^J'^^l^^, \/xieir, r,) K^^r) - ^^(r)] dr , (5-18) 

assuming that < /x < 1/2 and < u < 1. We returned to the integration over r in 
equation (5-15), and used the definition (5-5) of Q(t) and the expression (5-6) for a{r). We 
have also replaced uj^{r) by m^Q^(r), which is correct to lowest order in e(r) according to 
the definition (4-3) of the latter parameter. [Within the studied approximation, the same 
should be done with a;^(rj) involved in $j(rc, a); see the second of the equalities (5-14).] 

The integer n in the equation (5-16) should not be negative for the following reason. 
It is straightforward to see from equation (5-14) that in any case |Z_/Z+| < 1 (recall that 
< ^ < 7r/2). Hence 

|$(rc, a)\<- arctan (tan ^ | < - = — - < - . (5-19) 
' ^ ~ TT V 2 y - 2 4(1 -/x) 2 ^ ^ 

Obviously, /(rc, a) is non-negative for all Tc > r^, therefore the whole left hand side of 
equation (5-16) is always larger than —0.5 and thus cannot equal the right hand side with 
n = -l,-2,.... 

For a given set of parameters [including and j satisfying the selection rule (4-6)] and 
< a < 1 let us now study the roots of equation (5-16). First of all, e(r) = for a = 0; thus 
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I{rc, 0) = and no roots can be found. There is no c-mode for a strictly non-rotating disk. 
So the rotation parameter interval is < a < 1. Then it is easy to see from the expression 
(4-14) for e and the inequalities (2-2) that for any r > r^, 

e(r,re) < e^{r) = lim e(r,re) = ^ , ^ . (5-20) 

Therefore there exists a finite limit 

1 -y(r)r^+'^ / 

7(oo,a) = lim /(rc, a) = - / —J xi^oo(r) [m'^Q'i{r) - K'^(r)]dr 

VX, r 7(r)r^+- / [0(r)-0x(r)]K0i(r)-^^(r)] ^ 

Indeed, the integral (5-21) converges at the upper limit because its integrand is oc r^-'^'^/^ 
{m^ — 1) and oc r*^"^/^ (m^ > 1) when r — > oo [see the expressions (2-1)], and u is always 
less than one. For the special case a — 1, fi±(r) has a simple root at r = n = 1, so the 
singularity at the lower hmit becomes (r — rj)~('*+^/^). However, the integral still converges 
since // < 1/2. 

Hence I{rc, a) is a bounded function of Tc on the semiaxis, and there exists 

Im{a) = sup /(rc, a) , < /^(a) < +oo . (5-22) 

rc>ri 

By the inequality (5-19), we have already shown the existence of 

^m{a) = sup $(re, a) , $^(a) < 0.25(1 - /i)-^ . (5-23) 

rc'>ri 

This imphes that the left hand side of equation (5-16) is a bounded function of Vc'. 

Tn,{a) = sup J^(rc, a) , -0.25(1 - /x)"^ < J^rn{a) < +oo . (5-24) 

rc>ri 

But then equation (5-16) has no roots for all integers n > J-'m{o,), so only c-modes with the 
radial numbers < n < J-m{o,) can exist. In other words, for any value of a, < a < 1, 
only c-modes with the radial numbers n — 0,1,2, ... , N{a) can exist, with the integer N{a) 
satisfying J-m{o) — 1 < N{a) < !Fjn{o)- Along with the analytic dependence of J^{rc,a) on 
rc, this also shows that the set of roots of equation (5-16) for each admissible value of n is 
discrete; but we actually can do much better than that. 

Namely, we first recall from section 4.2 that for any positive value of a, the difference 
[Jl(r) — f^±(r)] is a decreasing positive function of r. Then the equalities (4-14), containing 
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the term — [fi(rc) — ^±{rc)], show us that both |a;(r, rc)| and e(r, Tc) are increasing functions 
of rc for any r > ri. This allows us to conclude that for any < a < 1 the left hand side 
J-{rc, a) of equation (5-16) is an increasing function of re- 

Indeed, from the equalities (5-14) and definitions (5-12) and (5-10) we see that for 9 ^ 
0, 7r/2, the phase $(rc, a) depends on Tc only through the combination |a;(rj, rc)|^ [^(n, ^c)]'', 
and the dependence is monotonically increasing. Since this combination is an increasing 
function of r^ so is the function $(rc, a) itself. In particular, its maximum is 

$r„(a) = $(oo, a) = lim $(a, rc) (5-25) 

[see equation (5-23)], whose value can easily be calculated explicitly according to equations 
(5-14), (5-12), (5-10) and (4-14). For the cases = 0, 7r/2, i. e., the boundary condition (5-1) 
of the second and first kind, respectively, $ is independent of Tc and a, 

$(rc, a) = ±— = constant , (5-26) 

4(1 ~ A*) 

and the relation (5-25) is still formally true. 
Next, the definition (5-18) shows that 



dl{rc,a) 1 1"'"= ^(r)r^+^ j xi[m'^^']_{r) — K,'^{r)] de{r,rc) , 



drc 27r i^. {r - ^ e{r,rc) drc 

so I{rc, a) is also an increasing function of Tc, with the maximum value [see (5-22)] 

/^(a) =/(cx),a) (5-27) 

given by expression (5-21). Thus ^(rc, a) = $(rc, a) -|- /(rc, a) increases with Tc to its 
maximum value at Tc = oo, 

J^m{a) = $m(a) + I,a{a) = $(oo, a) + /(oo, a) = J^(oo, a) . (5-28) 

Our last observation regarding J-{rc, a) is that 

-0.25(1 -/i)~^ < :F(ri, a) = <l>(ri, a) < 0.25(1 -/x)"^ < n , n = l,2,3,.... (5-29) 

The implication of the established properties of J-{rc, a) is clear and most important. 
Namely, for any given a (0 < a < 1), and any n — 1,2,3..., N{a), with the integer N{a) 
satisfying 

Tm{a) - 1 < N{a) < J^^(a) , (5-30) 
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there exists exactly one root rc (a) e (n, oo) of the increasing function [T{rc,a) —n], i. e., 
of equation (5-16) . Moreover, if 

J^(ri,a) = $(ri,a) < , J>„(a) = /(oo, a) + $(oo, a) > , (5-31) 

then in addition there exists a root r^\a) G (r^, oo) of equation (5-16) with n = 0. All the 
roots prove to be numbered in an increasing order, 

n < [rf (a)] < r^^\a) < rf (a) < . . . < rf \a) , N = N{a) . (5-32) 

Turning now to the eigenfrequencies as defined through the corresponding capture zone 
boundaries by the general formula (4-11), we can reformulate our result as follows. 

For any given a, < a < 1 the spectrum of the studied c-modes consists of a finite 
number of eigenfrequencies (T^"'^(a), n — [0], 1, 2, 3 ... , N(a); with N(a) specified by the in- 
equalities (5-30), and a^'^\a) existing only under the conditions (5-31). The eigenfrequencies 
are negative for m > and positive for m < 0, and their magnitudes are always ordered as 

\m\ [Qin) - Qx(n)] > [k(°)(a)l] > W^'\a)\ > k^'H^)! > ■ ■ ■ > W^'^Ka)] > , (5-33) 

with N = N{a). Note that this is true for a given value of the parameter F and any angular 
and vertical numbers m and j satisfying the selection rule (4-6), so that the complete notation 
for the eigenfrequency is (j^](a). 

Note also that if J-'m{0') < 1, then we have N{a) = by the right inequality in equation 
(5-30), and the only c-mode that can possibly exist is that with n — 0. On the other hand, 
the number N{a) is bounded for all values of a e [0, 1] by the maximum value 

^max = max J^m{a) = J^m{a*) , a* e [0, 1] , (5-34) 

0<a<l 

N{a) < Tjnax- So the total number of different c-modes with given angular and vertical 
numbers is less than T^ax + 1 for the whole family of corotating disks which differ only by 
the value of the angular momentum of the central body. If in particular Tmax < 1, which 
can happen if the speed of sound at the midplane is everywhere sufficiently large inside the 
disk, then at most one c-mode with n = may be excited in any disk of the family. 



5.4. Dependence of the spectrum on the rotation pctrameter: the cutoff values 

of a (n > 1) 

We now change our point of view: instead of studying all the roots for a given a, we 
specify some n>l and investigate how the trapping radius r^\a), determined by 

^ {ri-\ a) = $(rf , a) + Iiri^\ a) ^ n , (5-35) 



-18- 



behaves as a function of a. An intuitive expectation is that it monotonically decreases with 
a, but that might not be generally true. The differentiation of equation (5-35) gives 



da [dJ^/drc) 



=rc (a) 



and, although the denominator here is positive, it is difficult to establish that the sign of the 
numerator is also always positive. Even the sign of dl /da is not always obvious, because of 
many dependences on a in the expression (5-18): the quantities e, Oj^, and all depend 
on a. 

Nevertheless, /(cxd, a) goes to zero when a — > +0 (since ^ Q), so /(rc, a) < I{oo, a) 
goes to zero uniformly in > rj as well. In particular, /(r^, a) < n — 1/2 for small enough 
positive values of a. Since the first term on the left of equation (5-16) is less than 1/2 by 
equation (5-19), equation (5-16) has no root at all for such values of a. Therefore ri"'\a) 
does not exist in some positive vicinity of a = 0. 

Suppose now that ri"\a) docs exist for some a* > 0; that is, n < N{a^). By making 
the value of a smaller and smaller we would eventually find ri^\a) growing, to make up for 
the reduced value of I{rc^\a),a) < I{oo,a) +0. Recall that for a given a, I{rc,a) is a 
growing function of rc- At a certain value a„ > 0, only r^\an) — oo would satisfy equation 
(5-35), 

T (oo, ttn) — $(oo, a„) -|- /(oo, ttn) — n . (5-36) 
For a < a„, no value of the capture zone radius would satisfy it, so ri"\a) no longer exists. 

Therefore we have proved that for any admissible radial mode number n > 1 there 
exists a cutoff value a„ > such that ri^\a), and thus the corresponding c-mode, exists for 
a > an, and does not exist for < a < a„. When a — > a„ -f- 0, the capture zone boundary 
Tc (a) oo and the eigenfrequency — > 0. The cutoff value is determined by equation 
(5-36). Evidently, a„ increases with n, 

< ai < a2 < . . . . 

Also, the maximum radial mode number N{a) defined by the inequalities (5-30) is a piecewise 
constant increasing function with unit jumps at a — an- 

An approximate expression for the cutoff values can be obtained from equation (5-36) 
if they are small enough, so that a„/r"^^^(a„) <^ 1. We find that 



n + 0.25/(1 -/x) 



2 



, a^^ < sin 6' > ; 



n - 0.25/(1 -/x) 



2 



, ^ = ; (5-37) 
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where 

I-^l V-6)M V^^-l + -^^^ /(oc,a) = V^/+... . (5-38) 

We have used the fact that ri{a) ^ rj(0) = 6 for small a, along with the useful asymptotic 
expressions imphed by equations (4-14) and (2-1): 

6(r, re) = 2ar'/' (^^ - + . . . , ^(r, oo) = |^ + . . . ; (5-39) 

m^n\ - = m^ + - + ... (5-40) 



a;(r, Tc) = m 



1 ^ / 1 2 



c 



The latter, in turn, allow us to obtain, from equations (5-17), (5-14), (5-12), and (5-10), the 
results 

$(re,a) = -0.25(1 -/i)-^ + C» [a^''^^^-^'^) ,6^0; ^r^,a) = 0.25(1 -/i)-\ ^ = . (5-41) 

These relations, for = oo, arc also used to derive formulas (5-37). For the expressions 
(5-37) to be applicable, one essentially requires / ^ 1. 

The last item on our agenda for c-modes with n > 1 is to establish the behavior of 
rc^\a) and (^^"•'(a) near the cutoff value a„. Asymptotic calculations based on equations 
(5-35), (5-36), (4-14), and (4-13) allow us to find the results, which differ slightly depending 
on the angular number of the mode. 

a) Fundamental c-mode (m^ — j — 1). In this case we obtain, for a — > a„ -|- 0: 



r(")(a)«.C„(a-aJ-^, C„ ^ <> ^ll"? ^""^'^rZ ^^'^l (5-42) 

(7/4 - u) [$ (oo, an) + I (oo, an)\ J 

where prime denotes the derivative in a, 7oo is defined in section 5.1 [see equation (5-6)], and 
the expression (4-10) has been used. According to equation (4-13), for the eigenfrequency 
near its cutoff we have 

a(")(a) = ^ - (a - a^)^ ; (5-43) 



C3 

and the minus (plus) sign is taken for m = -|-1(— 1) 



b) C-modes with higher axial mode numbers (m? > 1/ Now we find (near the cutoff a 
a„ + 0): 



574= 



r(")(a)^D.(a-a„)-vi^, D^ = { ^^f?- ; (5-44) 

(5/4 - v) (00, an) + / (00, a„)J 



-20- 



with the value of xi given by equation (4-8). The corresponding eigenfrequency becomes 

2ma 2ma 



= - r f'r^'nS ~ (« - . (5-45) 



rr(a) 



n 



Recall from formula (5-26) that $'(oo, a) = for 9 = 7r/2, 0; that is, for boundary conditions 
of the first and second kind at the inner edge of the disk. Otherwise, the expressions for the 
coefficients C„ and Dn can also be made more explicit if a„ is small. 

Finally, we would like to say a few words about the behavior of the studied c-modes in 
the limit of very rapidly corotating disks, a — > 1 — 0. It turns out that the basic assumption of 
our approach becomes invalid in this limit, because the (assumed small) parameter e acquires 
large values at the inner edge of the disk, reaching eventually to infinity when a = 1 [unless 
the capture zone shrinks, i. e., unless ri^\a) rj(a) — 1 -|- for a — > 1 — 0]. That happens 
due to the fact that the denominator of the expression (4-3) or (4-14) for e goes to zero at the 
inner edge, ft±{ri = 1) = for a = 1. However, these misfortunes occur very close to a = 1, 
namely, for a > a±^ 0.953, at which value Jl±(r) as a function of r becomes non-monotonic 
by developing a maximum at some radius close to the innermost one. Therefore our formulas 
and conclusions remain valid for < a < a^, almost all the way up to the black hole limit 
a = 1. 



5.5. C-mode for n = 

Turning to the mode with the radial mode number n = 0, we rewrite the eigenvalue 
equation (5-16) as 

-$(re, a) = /(re, a) , n = . (5-46) 

Since I{rc,a) is positive, ^{rf..a) must be negative for the mode to exist. This condition, 
with the use of equations (5-17), (5-14), (5-12), and (5-10), translates into 

0.25xi(l -/x)-V(n)rf'^+'^e(r„r,)]'^ |a;(r,,r,)|' < ^tan^ . (5-47) 

Thus, the n = c-mode does not exist at all if = 0, i. e., for the boundary condition of 
the second kind at the inner edge of the disk. We consider now only positive values of 6. 
For such values, however, the mode is absent if the parameter a is small enough, because 
the right hand side of equation (5-46) is (9(a^/^) by equation (5-38), while the left hand side 
approaches 0.25(1 — fx)^^ by equation (5-41). A sufficient condition for the mode to exist at 
least for some values of a is evidently 

max I(oo, a) > —j- r . 

0<a<0.95 ^ ' 4(1 — 11) 
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If it exists, the n — mode thus has a cutoff value ao, as the modes with other radial 
numbers do, and this cutoff is the smallest one. Its approximate magnitude is given by the 
first of equations (5-37) with n — 0, 

1 



16(1 -n^P ' 

and the behavior of the capture zone boundary and the eigenfrequency near the cutoff is 
described by either expressions (5-42), (5-43) or (5-44), (5-45), with n set to zero. 

Therefore there arc no c-modes of the studied type when the disk rotation is sufficiently 
slow. Namely, if the boundary conditions at the inner edge are of the second kind, then there 
are no c-modes for < a < oi. For all other boundary conditions, the same is true in a 
smaller interval of rotation parameter values, < a < Oq < Oi . 



This concludes our theory of c-mode oscillations of accretion disks. Its quantitative 
accuracy is limited by the WKB approximation used throughout, as in previous papers 
(Nowak & Wagoner 1992; Perez et al. 1997). On the other hand, its qualitative results, that 
is, the above properties of the c-mode spectrum and its dependence on the rotation parameter 
a, rest on a single assumption that the speed of sound does not drop too fast with radius, 
which seems to be valid for existing models of accretion disks. If, nevertheless, such a fast 
drop does take place [i/ > 5/4 in our terms, see formula (5-6)], then the qualitative picture 
changes drastically. The set of the radial mode numbers becomes infinite, there are no cutoff 
values of the rotation parameter a, and the c-modes of the studied type are present in any 
corotating disk, no matter how slow its rotation. As it turns out, for i/ < 5/4 the same effect 
on the results is caused by using a very rough approximation of the radial eigenf unction, 
namely, just the Airy function (5-8) throughout the whole trapping zone (Silbergleit & 
Wagoner 1998): in this approximation the number of the radial modes becomes infinite, and 
the cutoff values of parameter a disappear, so that all the modes exist for any corotating 
disk. 



6. Numerical Results and Discussion 



We now apply the above analysis to models of black hole accretion disks. The numerical 
results discussed below were obtained by assuming that gas pressure dominates within the 
disk, so that F = 5/3. However, this is only true near the inner edge [justifying our choice 
H = 2/5 in equation (5-6) consistent with the results of Page & Thorne (1974)], at radii 
r ^ ri, and at all radii for luminosities L <^ LEdd- 
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Recall that the interior structure of the (zero buoyancy) accretion disk enters our formu- 
lation only through the constant F and the function a{r) (not to be confused with the usual 
viscosity parameter, here denoted by «*)• inversely proportional to the speed of sound on the 
midplane and modeled by equation (5-6). In this equation we also chose 7(r) to be constant, 
and took z/ = unless otherwise indicated. In section 4.3, we found that the properties of 
the modes do not change greatly when small amounts of buoyancy are introduced. 

For all of the tables (but not the figures), we used M = IO^Mq, = 0.01, and 
L = O.lLEdd so we could compare our results with those of Perez (1993). For a standard 
thin accretion disk (Shakura & Sunyacv 1973), at large radii (where the disk is nonrelativistic 
and gas pressure dominates) the speed of sound is given by 

corresponding to = 9/20 (Nowak 1992). (We are here reverting to ordinary units.) Since 
the physical conditions within the disk are more uncertain near its inner edge (but where 
gas pressure should also dominate), we have also used the nonrelativistic expression 

cjc = 4.32 X 10-='^^^^^[c^(r - r.)/GM]2/^ (r - u) (6-2) 

there (Nowak 1992), as indicated previously. We employ this relation (6-2) when using 
1/ — in equation (5-6). We obtain the (constant) value of 7 used in equation (5-6) by 
setting a{r) — l/cs{r) there (again neglecting the relativistic corrections near the inner 
edge). The various relativistic corrections only become significant when a > 0.5. We have 
only studied the fundamental c-mode (m^ = j = 1), but have employed various values of 
the radial mode number n. 

In Table 1 we present the cutoff values of a [an, from equation (5-36); for a < an, the 
mode does not exist] for radial mode numbers n = 0,1,2 and two choices of the boundary 
condition parameter. There is no c-mode for n = and = 0; no solution is also indicated 
by no entry in the subsequent tables. The results for 9 = 7r/4 are essentially identical to 
those shown for 9 — 7r/2. The results agree with the analytic expressions (5-37) vahd for 
values of a„ ^ 1. These results ior 9 — 7r/2 are extended to larger values of n in Table 2. 
The maximum n — 175 included corresponds to the requirement a < a± ^ 0.953 for the 
validity of our analysis. 

Table 3 presents values of our basic expansion parameter e{ri, Tc) , obtained from equation 
(4-14), as a function of a for the lower radial mode numbers and two choices of boundary 
condition. We see that it is indeed much less than unity, as also required for our analysis. 
Recall that e(ri, Tc) > e(r, Tc) > for all Vc > r > rj. The same format is used for Tables 4 
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and 5. The eigenfrequency a (in units of /GM, as usual) is presented in Table 4. We see 
that it is a decreasing function of the radial mode number n, more strongly for lower values 
of a. The corresponding radial extent r,, — r, of the mode (in units of GM/c^, as usual) 
is shown in Table 5. It increases significantly with radial mode number. We note that the 
sensitivity of the results in Tables 3-5 to the boundary condition parameter 6 is modest. 

Table 6 shows the dependence of the outer radius Tc of the capture zone on values of 
a near its cutoff a„, for n = and 9 = 7r/2 (corresponding to = 3.22 x 10^^). From 
these data, one can see that when log(a — Oq) ^ "~6, the asymptotic formula (5-42) holds: 
Tc oc (a — ao)~'^. The exponent q = 0.56 is very close to that predicted, q = 4/7, for u = 0. 

The effect of changing the behavior of the speed of sound at large radii is illustrated in 
Table 7. The radial size of the mode for u = (from Table 5) is compared to that obtained 
with the choice u = 9/20 of the standard accretion disk model, equation (6-1). This shows 
that the decrease in the speed of sound with radius also reduces the size of the trapping 
region, as expected. Correspondingly, the values of a„ shown in Table 1 are reduced by 
about a factor of 7.5, and the values of e(ri, Tc) in Table 3 are reduced by about a factor of 
2. The eigenfrequencies in Table 4 are increased by factors of 1.2-2.8 (with n = 0, 1) for 
a — 10~^, but were unaffected for a — 0.5, as expected. 

We present our major observationally relevant results in the following figures. For 
them, we have chosen the radial mode number n — and the boundary condition parameter 
9 — 7r/2. (One might expect that the lowest radial mode could be the one most easily excited 
and with the largest net modulation.) The speed of sound parameter is taken to be = 
for Figures 2 and 3, and u — 9/20 for Figures 4(a,b) (which differed very little when u — 
was used). 

In Figure 2, we plot the relation between the fundamental c-mode and g-mode eigen- 
frequencies (scaled by the black hole mass) and the black hole angular momentum. This 
illustrates the dramatically different dependence of the c-mode frequency. Also shown is 
the orbital frequency flmax = ^{f^i) of a (commonly invoked) 'blob' at the inner disk radius. 
As expected, the c-mode frequency approaches Qmaxf^T^ as a ^ 1. The case u = 9/20 lies 
almost entirely within the band containing the range of masses and luminosities indicated. 
It is significant that the c-mode results of Perez (1993), obtained by numerical integration of 
equations (2-4) and (2-5), agree within the band shown in this figure (obtained via a further 
radial WKB approximation) . 

In Figure 3, we present our numerical results (points) for the size of the trapping region. 
The curves shown are the fits of these results to the formula 



re - n = Xo(GM/c')a-^i(l - a) 



(6-3) 
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giving 

Kq = 0.058 , Ki = 0.66 , K2 = 0.31 (M = lOM©) , 
Ko = 0.021 , Ki = 0.55 , K2 = 0.39 (M = lO^M©) . 

This illustrates the fact that the radial extent of the mode is only appreciable for slowly 
rotating black holes (a <^ 1). To obtain the corresponding cigenfrequency, one can then use 
equation (6-3) and the known function = Mf{a) to obtain the value of the trapping radius 
Tc- Our fundamental result, equation (4-11) [or equation (4-13) for a -C 1] with m = — 1, 
then gives \a\. 

The c-mode frequency (scaled by mass) depends only on the black hole angular mo- 
mentum and the accretion disk speed of sound. The first dependence has been shown in 
Figure 2. Rather than showing the second dependence directly, it is more relevant to instead 
use an observable, the luminosity. For fixed a* and M, we use equation (6-1) to relate 
it to the speed of sound. We want to know how the frequency changes as the luminosity 
(proportional to the mass accretion rate) varies with time. That dependence is known to be 
very weak for the g-modes (Perez et al. 1997). In contrast, the m — fundamental p-mode 
has a frequency oc c^^ and a radial size (r_ - Ti) oc cj^ (Kato & Pukue 1980; Nowak & 
Wagoner 1991; Perez 1993). 

We fit the dependences of the eigenfrequency on accretion disk luminosity shown by the 
points in Figures 4(a) and 4(b) by the form \a\ oc L~^(^'"). The results shown correspond 
to 

K{10Mq, 10'^) = 0.30 , ^(lO^Mo, 10-2) = 0.039 ; 
K{10Mq, IQ-^) = 0.020 , K{10^Mq, 10"^) = 0.0044 . 

Note that the dependence is weaker for the larger value of a, again as expected since the 
properties of the disk (except its inner radius) become irrelevant as a — > 1. 

Although the fundamental c-mode is almost incompressible, the changing projected area 
of the mode could modulate the luminosity via reflection of radiation from the postulated 
'corona' surrounding the accretion disk. Of course, this requires that the disk not be viewed 
close to face-on. The observability of a c-mode induced modulation of the detected flux 
would seem to require large values of Vc — r^. This in turn would imply small values of 
a, from Figure 3, and correspondingly small values of frequency, from Figure 2. We then 
see from Figure 4(a) that the dependence of the frequency on luminosity is relatively weak 
but might be detectable for the stellar mass black holes. Issues such as the excitation and 
damping of the c-modes, including their leakage into the black hole via accretion from the 
inner edge of the disk, are beyond the scope of this paper. 
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Finally, we emphasize that of the fundamental (g, p, c) modes, only the c-mode can 
(generically) have a frequency \(t\ <^ Q.. The c-modes are candidates for those low frequency 
features in the power spectra of accreting black holes whose frequency varies only weakly 
with changes in luminosity. An example of a candidate is the 9 Hz modulation in the 
'microquasar' GRO J1655-40 observed by the RXTE satellite (Remillard et al. 1999). Since 
the mass of the presumed black hole in this binary has been determined to be M = IMq 
(Shahbaz et al. 1999), this frequency requires a black hole angular momentum a = 0.18 
(Figure 2) if it is produced by a low n c-mode. 

However, wc sec from Figure 3 that the lowest radial mode would have a radial extent 
of only Q.2GM / . (But from the results in Table 5, we expect the radial extent of the n = 1 
mode to be about 3 times greater.) In addition, the energy spectrum of the X-rays when this 
modulation was present was 'softer' (more thermal) than typical. This may complicate the 
above proposal that the photons are modulated via disk reflection from the corona. Clearly, 
one should search for other radial (or vertical) modes to confirm any identification. 

This work was supported by NASA grant NAS 8-39225 to Gravity Probe B and NASA 
grant NAG 5-3102 to RVW, who also thanks the Aspen Center for Physics for support during 
a 1999 summer workshop. We arc grateful to Lev Kapitanski for finding the number theory 
results for the Diophantine equation (4-6), and to Dana Lehr for her remarks and help with 
the figures. 
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Table 1. Cutoff values (a^) of a 
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Table 3. Expansion parameter e(rj,rc) 
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Table 4. Eigenfrequencies —a/{c'/GM) 
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Table 5. Trapping region [vc — ri)/{GM/c'^) 
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Table 6. Capture radius for a near its lower limit (n = 0) 
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Table 7. Dependence of - ri)/{GM/c^) on v 
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Fig. 1. — Radii between which o;^ < k^, as a function of the eigenfrequcney. Three choices 
of axial mode number m and black hole (dimensionless) angular momentum a are included. 
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Fig. 2. — The dependence of three characteristic disk frequencies on the angular momentum 
of the black hole. The top one is the orbital frequency of a 'blob' at the inner edge of the 
accretion disk. The c-mode curves bound the region corresponding to the parameter ranges 
0.001 < L/LEdd < 0.1 and 10 < M/M© < 10^ with a* = 0.1. 
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Fig. 3. — The dependence of the radial extent of the mode trapping region on the angular 
momentum of the black hole, choosing L = O.lLEdd- For M = IOMq (crosses), a^, = 0.1, 
while for M = IO^Mq (hexagons), a* = 0.01. The fits are of the form Koa~^^{l — a)^^. 
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Fig. 4. — The dependence of the eigcnfrequency on the accretion disk luminosity, choosing 
a* = 0.1. For plot (a), a = 0.001 and for plot (b), a = 0.1. 



